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ABSTRACT 

Observations of the black hole X-ray binary V404 Cyg with the very long baseline in- 
terferometer HSA (the High Sensitivity Array) have detected the source at a frequency 
of 8.4 GHz, providing a source position accurate to 0.3 mas relative to the calibrator 
source. The observations put an upper limit of 1.3 mas on the source size (5.2 AU at 
4kpc) and a lower limit of 7 x 10 6 K on its brightness temperature during the normal 
quiescent state, implying that the radio emission must be non-thermal, most probably 
synchrotron radiation, possibly from a jet. The radio lightcurves show a short flare, 
with a rise time of ~ 30 min, confirming that the source remains active in the quiescent 
state. 

Key words: X-rays: binaries - radio continuum: stars - stars: individual (V404 Cyg) 
- ISM: jets and outflows - astrometry 



1 INTRODUCTION in two syste ms, GRS 1915+105 (|Dhawan et alj |2000T ) and 

. , , , „ .. , . Cygnus X-l dStirling et al1l200ll) . at X-ray luminosities of 5 

Accreting black hole X-ray binary systems spend the ma- , „ „ rj.i_r-.jj- i i 

° , i, • anc ' U.o per cent ol the Eddington luminosity respectively, 

jority ol their time in the so-called quiescent state, 



a low-luminosity state (below a few x 10 of the Ed- 



With only two relatively luminous examples of self-absorbed 

jets, it has not been possible to study the dependence of the 
dington luminosity) whose X-ray emission is characterised . , . , , , 

° , _ ' . jet size and geometry with system parameters such as or- 

by a hard (1.5 < T < 2.1) power-law spectrum , ., , . , „ ,, , ,. , 

i j, „, ; — r— — n I 1. ' , , , bital period, nature ol the mass donor, accretion rate, and 

(|McClintock fc Remillardl 12009 .. Despite such highly sub- , , . , ' 

black hole mass. 

Eddmgton systems being so numerous, their intrinsic tamt- 

ness poses serious observational challenges to any attempt to On moving from the hard to the quiescent state, the ra- 

understand the nature of the accretion flow in these systems. dio properties of X-ray binaries become progressively less 
Accreting black holes have been studied in much more constrained, owing to the limitations of current instru- 
depth in the "hard" state, a regime with a similarly hard mentation. In particular, the questions of whether radio- 
spectrum but occurring at higher X-ray luminosities (up to a emitting outflows persist into quiescence, and if so, what 
few per cent of the Eddingt on luminosity ). In this state, they form they take , remain to be answered. Quie scence is cur- 
show flat radio spectra, (e.g I Fen der 200l]) interpreted as syn- rently defined (|McClintock fc Remillardl 120061 ) only by the 
chrotron emission arising from part i ally self-absorbed, con- luminosity criterion L x < 10 33 ' 5 ergs _1 (i.e. Lx/iEdd < 
ical j ets (jBlandford fc Konigll fl979l : iHiellming fc Johnston! 2.3 x lO~ e (A/ B H/lOAf ), where L E dd is the Eddington lu- 
1 19881 ). Such jet-like structures have been directly imaged minosity and M B h is the black hole mass), although there 

are indications that the X-ray spectrum in quiescence soft- 
ens from the canonical powe r-law index F ~ 1 . 5 of the hard 
* email: jmiller@nrao.edu state to a value closer to 2 (|Gallo et alj 12009 : ICorbel et alj 
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l200d : iBradlev et al.|[2007h . However, at such low luminosi- 
ties, instrumental sensitivity begins to hamper efforts to 
characterise the behaviour of quiescent systems, and the 
detailed spectral properties of the quiescent state remain 
to be determined. Extrapolation of the observed hard state 
behaviour suggests that it is a reasonable assumption that 
jets should also exist in the quiescent state. This hypothesis 
is supported by the detection of A0620- 00 in the quiescen t 
state in both the radio and X-ray bands (|Gallo et al. 
which demonst rated that the rad io/X-ray correlation seen in 
the hard state (|Gallo et alj|2003h remains unbroken down to 
10 -8,5 Z/Edd- Further support is provided by the o bserved flat 
radio spectrum in V404 Cyg (|Gallo et al.| [2005). However, 
the evidence is circumstantial, and direct empirical proof 
for the existence of such quiescent jets is still lacking. 

1.1 V404 Cyg 

V404 Cyg (GS 2023+338) is the most luminous of the qui- 
escent black hole X-ray binaries, with a mean X-ray lu- 
minosity of ~ lO^friAlkpc^ergs" 1 dGarcia et"all l200ll ; 
iKong et al.ll2002l ; iHvnes et aj]|2004l ; IBradlev et al.ll2007t l. Its 
high, accurately-determined mass function, 6.08 ± O.O6A/0 
(the most statistical ly significant measureme nt of a mass 
function above 3 Mm ; ICasares fc Charledfl994l ), implies that 
the compact object must be a black hole. The best-fitting 
syste m parameters indic ate that the mass donor is a K0 IV 
star (|Casares et al.lll993|) in a 6.5-d or bit with a black hole of 
mass 12 ± 2Mq ( Shahbaz et al.|[l994l '). although contamina- 
tion of the infrared spectrum by light from the accretion disc 
could reduce this to I O Mq , with a firm lower lim it of 7AMq 
jShahbaz et al.lll996l ). lHan fc Hiellminsl (| 19921 ) measured a 
kinematic Hi distance of d ^ 3.2 kpc. Assuming Eddington- 
limited X-ray emission during the 1989 outburst gives a 
maximum distance of ^ 3.7 kpc, whereas the optical mag- 
nitude and reddening info rmation suggests d = 4.0l \% kpc 



ijjonker fc NelemansH20o3 ). We will adopt the latter value 
for the rest of this work. 

The radio spectrum of V404 Cyg was measured to 
be flat, with a spectral index of a — 0.09 ± 0.19 (Gallo, 
Fender & Hynes 2005), defining S v oc u a . Its flux density 
is ~ 0.3 mjy at radio wavelengths (Gallo, Fender & Pooley 
2003), although it is k nown to vary up to levels of 1.5 mjy 
jHiellming et alJl200Ch . The radio emission is believed to be 
of synchrotron origin ijGallo et alJl2005l ). and the quiescent 
radio and X-ray luminosities are consistent with an extrap- 
olation of the non-linear r adio/X-ray correla tion found for 
hard state X-ray binaries l)Gallo et al.l 120031 ). suggestive of 
jet emission. However, the collimated nature of the outflow 
rem ains to be proven. 

iMuno fc Mauerhanl (|2006l l detected 4.5- and 8- /im in- 
frared emission from V404 Cyg in excess of that expected 
from the Rayleigh- Jeans tail of the companion star, which 
they attributed to emis sion from the accretion disc. How- 
ever, iGallo et al] (|2007h demonstrated that the broadband 
radio through X-ray spectral energy distribution (SED) 
could be well fitted with a maximally jet-dominated model, 
suggesting that the excess infrared emission was consis- 
tent with being synchrotron emission from a partially self- 
absorbed jet. 

As the most luminous quiescent black hole X-ray binary, 
and with so much circumstantial evidence for jets, V404 Cyg 



is the ideal candidate in which to try to demonstrate empir- 
ically the existence of jets in a quiescent accreting black hole 
system. Since the source is so faint (~ 0.3 mjy; other qui- 
escent systems are fainter still), we proposed for Very Long 
Baseline Interferometry (VLBI) observations with the High 
Sensitivity Array (HSA) in an attempt to resolve the jets, 
and as a pilot to check the feasibility of a future parallax 
measurement to determine the source distance. 



2 OBSERVATIONS AND DATA REDUCTION 

V404 Cyg was observed for 4.25 h with the HSA on 2007 De- 
cember 2, from 18:00:00 until 22:15:00 UT, under program 
code BG 168. We used nine antennas of the Very Long Base- 
line Array (VLB A; the tenth antenna, St. Croix, was down 
for maintenance) plus the Green Bank Telescope (GBT), Ef- 
felsberg and the phased Very Large Array (VLA). Eleven of 
the 24 available VLA stations were retrofitted EVLA anten- 
nas. We observed at 8.4 GHz, in dual-circular polarization, 
with a bandwidth of 32 MHz per polarization. Samples were 
quantized to 2-bits, giving a total bit rate of 256 Mb s _1 . We 
used J 2025+3343, a 2-Jy source located only 16.6 arcmin 
from the target, as both the fringe finder and phase refer- 
ence source. We observed in 3-min cycles, spending 2min 
on the target and 1 min on the calibrator in each cycle. Ev- 
ery 24min, we also made observations of a check source, 
J 2023+3153, the next nearest bright calibrator (1.9° from 
the phase reference source), followed by a rephasing of the 
VLA. 

Data were reduced according to standard procedures 
within the National Radio Astronomy Observatory' s Astro- 
nomi cal Image Processing System (AIPS) software (jGreisenl 
2003). During the self-calibration of the phase reference 
source, the initial model of the calibrator was created using 
only the VLB A antennas, since the a priori calibration (us- 
ing system temperatures) for the VLBA antennas is better 
than that for the larger dishes. Amplitude self-calibration 
was then used to correct the gains of the large dishes ac- 
cording to the VLBA-only model. Since the sensitivites on 
the typically long baselines connecting the large dishes dom- 
inated the self-calibration solutions, their absolute flux scale 
could not be sufficiently well determined using standard cal- 
ibration. The visibilities on baselines from the Los Alamos 
(LA) station to the VLA were compared to those at simi- 
lar (u,v) co-ordinates between the LA and Pie Town (PT) 
stations, and a global scaling factor of 0.89 was applied to 
correct the VLA gains and bring the two sets of visibilities 
into agreement. A similar procedure was applied to baselines 
from the North Liberty (NL) station to both the Hancock 
(HN) antenna and the GBT, although it was found that no 
scaling was necessary for the GBT gains. Finally, the calibra- 
tion derived on the phase reference source was transferred 
to the target. Owing to the low flux density of the target, 
no self-calibration was performed. 

Th e extreme scat tering towards the phase reference 
source (|Fev et al.lll989l ) meant that the achievable resolution 
was set by the longest unscattered baseline to the calibrator, 
rather than the expected instrumental resolution. For this 
reason, the data from Effelsberg had to be rejected. Visibil- 
ities were scattered to zero amplitudes on baselines longer 
than 150 MA (Fig. [TJ, whereas baselines to Effelsberg were 
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Figure 1. Amplitude of the real part of the visibility plotted 
against radial distance from the centre of the ui>-plane, for the 
phase reference source J 2025+3343. The shape of the plot is con- 
sistent with anisotropic int erstellar scattering (as parameterised 
by, e.g., |Pesai &: Fe"y||200lf ) giving rise to an elliptical scattering 
disc of size 2.2 X 1.4 mas in PA 2° (see Section[5}, with the power- 
law index of the electron density fluctuations in the i nterstellar 
mediu m being close to 4, as found for J 2025+3343 bv lFev et al.l 
dl989h . 



between 120 and 290 MA. Those in the range 120-150 MA 
were taken at the lowest elevations, and it was not deemed 
possible to achieve accurate gain calibration for any of the 
data from this station. 

The VLA data were independently reduced and im- 
aged, using standard procedures within AIPS. Weather con- 
ditions at the VLA were excellent throughout the observa- 
tions. Amplitude calibration was performed using observa- 
tions of 3C48, adopting the flux density scale set by the 
coefficients derived at the VLA in 1999.2, as implemented 
in the 31DEC08 version of AIPS. While no proper polar- 
ization calibration could be performed, the low polarization 
leakage of the VLA (of order a few per cent) meant that we 
could obtain approximate constraints on the linear (P) and 
circular (V) polarization of the source by imaging the two 
internal IF pairs separately in all four Stokes parameters I, 
Q, U and V, and forming P = (Q 2 + U 2 ) 1/2 . 



3 RESULTS 

V404 Cyg was detected by the HSA at a 12a level with a 
mean flux density of 0.32 ± 0.03mJybm _1 (here and in all 
the following sections, unless specified, we quote la error 
bars). The source was unresolved down to a beamsize of 
5.0 x 1.4 mas 2 in PA 2?0 (Fig. [2j. The elongation is caused 
by the high sensitivities of the GBT and phased VLA, with 
baselines to (and particularly between) those antennas dom- 
inating the naturally-weighted sampling function. 

We independently reduced the VLA data, to measure 
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Figure 2. Image of V404 Cyg. The beamsize (5. Ox 1.4 mas 2 in PA 
2?0) is shown in the lower right. Contours are at levels of ±(v / 2)" 
times the rms noise in the image (0.027 mjy bm -1 ), where n = 
3,4,5,..., with the lowest contour being at 0.078 mjybm -1 , and 
the peak flux density in the image being 0.322 mjy bm — 1 . The 
source is unresolved. 



the total flux density of the source and ascertain whether 
any emission was being resolved out by the HSA. The flux 
density measured by the VLA, averaged over the length of 
the observation, was 0.47 ± 0.03 mjy, showing that the HSA 
is detecting most of the flux density, and constraining the 
majority of the emission to come from a region smaller than 
5 mas (Section|4]discusses the size constraint in more detail). 



3.1 Variability 

To search for shorter-timescale variability, the VLA data 
were split into short (30-min) time segments and im- 
aged. The source was fit in the image plane with a point 
source model during each 30-min interval, and the resulting 
lightcurve is shown in Fig. [3] 

The measured HSA flux density is consistent with recov- 
ering ~ 70 per cent of the total flux density measured with 
the VLA at all times. When imaging the data from the flare 
at the end of the observations, editing out the data from the 
phased VLA resulted in recovering a much higher fraction of 
the total source flux (1.00 + 0.14 mjy out of 1.14+0.07 mjy), 
suggesting possible amplitude calibration or phasing prob- 
lems with the VLA. However, a similar analysis performed 
on the non-flare data showed no significant discrepancy be- 
tween the results with and without including the phased 
VLA. Thus any problems with the phased VLA were con- 
fined to the time of the flaring observations. 



3.1.1 Flaring emission 

The source flux density as measured by the VLA rose 
by a factor of 3 during the flare at the end of the ob- 
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Figure 3. VLA (filled circles) and HSA (open triangles) 
lightcurvcs during the HSA observations of V404 Cyg. Note the 
flare at the end of the observations. 



serving run, and there is also the suggestion of a smaller 
flare earlier in the observation, at 18:45. The rise time 
of the main flare was < 30min, and it showed no de- 
tectable circular or linear polarization (to a 5a limit of 
32 per cent), ruling out coherent stellar bursts from the 
donor as the origin of the emission. However, incoherent 
gyrosynchrotron stellar flares with timescales of tens of 
minutes and much lower degrees of circular polar ization 
are c ommon in RS CVn and Algol-type systems (|Giidell 
120021 ), and the peak radio luminosities of such flares can 
reach levels of 10 17 -10 18 ergs" 1 Hz" 1 (e.g.lMutel et al.lll987l ; 
iBenz fc Gude3ll994l ; iGarci'a-Sanchez et al.ll2003l ). One such 



system, UX Ari, with a K0 subgia nt in a 6.44-d orbital pe - 
riod with a G5 main sequence star (jStrassmeier et aTll 19931 ). 
makes for a good comparison to V404 Cyg in terms of the 
rotation rate and stellar type of the flaring subgiant star. 
The maximum flare strength seen in this sy stem in 3 years 
of ra dio monitoring was 0.82 Jy at 8.3 GHz (Richards et al. 
120031 ). implying a peak luminosity of 2.5 x 10 18 erg s 1 Hz 1 , 
among the highest observed to date in such systems. The 
peak luminosity of the flare seen in V404 Cyg was an order 
of magnitude higher still, making it unlikely (if not impossi- 
ble) that the observed flare in our data originates from the 
donor star. 

Splitting the HSA observations into 1-h chunks for 
imaging (Fig. |3}, we found that the flare was also seen in 
the high-resolution data. Imaging only the data from during 
the flare also showed that the source was unresolved, being 
smaller than the beamsize of 4.7 x 1.1 mas 2 in PA — 5?9. 



3.1.2 N on- flaring emission 

Taking only the data prior to the flare (18:00-21:30 UT), 
the HSA measured a mean quiescent flux density of 0.24 ± 
0.03 mjy beam" 1 , and the source was unresolved with a 
beamsize of 5.2 x 1.4 mas 2 in PA 5?5. The VLA alone mea- 
sured 0.35 ± 0.03 mjy beam" 1 . There was no evidence for 
any source polarization with the VLA, down to a 5a limit 
of 40 per cent for both linear and circular polarization. 



3.2 Source position 

The shift in source position between quiescent and flar- 
ing states was 0.13 ± 0.31 mas (la uncertainty), consis- 
tent with zero. Using all the available data, the mea- 
sured VLBI position, relative to the phase reference source 
J 2025+3343 (whose position was taken to b<£] (J 2000) 
20 h 25 m 10!8421050(224) 33°43'00'.'21443(42)) was 

RA = 20 h 24 m 03 l :821269 ± 0.000004 

Dec. = +33°52'01"8988 ± 0.0002 (J 2000), 

where the quoted error bars are the purely statistical errors 
from fitting in the image plane. Imaging of the check source, 
J 2023+3153, phase referenced to J 2025+3343, showed that 
the shift in its position, relative to the catalogued posi- 
tion, was 0.34 mas in RA and 0.09 mas in Dec, compared 
to uncertainties of 0.24 and 0.41 mas respectively in the cat- 
alogued position of J 2023+3153. Since the check source was 
seven times further away from the phase reference calibra- 
tor than the target source, we are therefore confident that 
the systematic errors in the source position (e.g. due to at- 
mospheric phase gradients) are small and the positional un- 
certainty is dominated by the signal-to-noise ratio on the 
target. 

The unresolved nature of the source makes this system 
a good candidate for a future parallactic distance determina- 
tion, free from the extra positional uncertainties introduced 
by source structure. The current measurement of the posi- 
tion would b e used as an initial startin g point. A distance of 
4.0l 2 ;§kpc (jjonker fc Nelemand I200II ) would imply a par- 
allax Of 7T = 0.2510.U8 mas, which is certainly feasible with 
current instrumentation. 



4 JET SIZE CONSTRAINTS 

4.1 Image and visibility- based limits 

In order to put limits on the size of any extension, we fit- 
ted the source with an elliptical Gaussian in both the image 
plane (with the AIPS task JMFIT) and the u^plane (us- 
ing the AIPS task UVFIT and the Caltech VLBI package 
Difmap). Using the whole data set with all available sta- 
tions, the (la) upper limit on the FWHM of the major axis 
was found to be ~ 1.3 mas. We note that a sufficiently low- 
level extension would not show up in the fitting results, so 
this is not a rigorous upper limit to the source size. 



4.2 Brightness temperature limits 

Radio surface brightness may be parameterized in terms of 
the temperature a blackbody would need to have in order 
to produce a source of the observed size and flux density, at 
the observed frequency, in the Rayleigh- Jeans limit. This is 
the brightness temperature, given by 

where c is the speed of light, _?„ the source flux density at 
frequency v, A?b is Boltzmann's constant and the solid an- 
gle subtended by the source. Our most stringent constraint 

http: / /vlbi. gsfc.nasa.gov/solutions/2005f_astro/2005f_astro_cat.txt 
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comes from the microflare, when the flux density measured 
by the HSA was 1.0 ±0.1 mjy and fitting an elliptical Gaus- 
sian to the source in the image plane gave an upper limit to 
the source size of 1.8 x 1.7 mas 2 . This gives a lower limit to 
the brightness temperature of Tb > 1.1 x 10 7 K. For compari- 
son, for the non-flaring data, fitting an elliptical Gaussian in 
the image plane gave an upper limit to the size of 2 x 0.6 mas 2 
and a consequent lower limit to the brightness temperature 
of 7 x 10 6 K. 

Previous VLBI observations of Galactic X-ray binary 
jets have measured brightness temperat ures of ^ 10 9 K (e.g 
iDhawan et"al]|2000l : lMirabel et al.ll200lh . Should this be the 
case in V404 Cyg, we can use the flux density during the 
flare to place an upper limit of 4.6 x 10 -19 sr on the solid 
angle subtended by the source, corresponding to a geomet- 
ric mean size of 0.16 mas. However, most resolved X-ray bi- 
nary systems have shown collimated jets, so as an example, 
for an axial ratio of 4:1 (the definition of jets adopted by 
IBridle fc Perieylll984l ), the derived solid angle would trans- 
late into a jet length of 0.3 mas. 

The inverse Compton catastrophe limits the bright- 
ness temperature of a radio-emitting source to < 10 12 K, 
if it is not relativistically boosted. For the flux density of 
1.0 mjy beam - measured during the flare, this limits the 
source size to > 5 /^as (> 2 /xas for the non- flaring emission 
with a flux density of 0.24 mjy beam -1 ). 

4.3 Variability limits 

The < 30-min rise time of the flare corresponds to a distance 
of 5.4 x 10 1 m for a source moving at the speed of light, 
or an angular size of 0.9 mas if the motion is in the plane 
of the sky. However, sources moving at small angles to the 
line of sight can show apparent superluminal motion, with 
/?app > 1. We have no constraints on the speed or inclination 
angle to the line of sight of the source, but constrain the size 
of the flaring region to be < 0.9/3 app mas. 

Similarly, for the brightness temperature size con- 
straints derived in Section [4.21 we can derive limits on the 
jet speed, for a rise time of 30min. For Tb ^ 10 9 K, the jet 
speed is limited to < 0.33c, and for Tb < 10 12 K, the jet 
speed is > 830 kms -1 . If the jet is to be moving close to the 
speed of light, it would have to be more elongated, and have 
a lower brightness temperature. 

4.4 Theoretical limits 

iHeinzl (|2006T ) derived an expression for the length of a fiat- 
spectrum, partially self-absorbed jet, as a function of basic 
jet parameters, normalizing the relati on using the VLBA 
measurements of the jet in Cygnus X-l (jStirling et alj|200ll) . 
The observed angular extent of the jet, /x 50 %, within which 
50 per cent of the flux density is contained, can be written 
as 

^g(sine) 10 V /17 
1*50% oc I ^ J mas, (2) 

where d is the source distance, 8 — [r(l — /3cos(9)] -1 is the 
Doppler factor, 9 the inclination angle of the jet to the line 
of sight, T the Lorentz factor (1 — /9 2 ) -1//2 , and P = v/c is 
the jet velocity in units of the speed of light. IShahbaz et al.l 



1 19941 ) derived an inclination angle of 56 ± 4° for V404 Cyg , 
and we measured a flux density of 0.3 mjy. IHeinzl (|2006l ) 
took values of (J,so% = 3 mas, S v = 12 mjy and 6 = 35° 
for Cygnus X-l. With these parameters, a distance of 4kpc, 
and assuming that the Doppler factor is similar in the two 
sources (since /iso% scales as S~ 2 ^ 17 , the exact value is unim- 
portant) and that the disc is perpendicular to the jet, we 
predict a jet length of 0.6 mas in V404 Cyg. 

4.5 Previous constraints 

iMioduszewski et all (|2008l ) also observed V404 Cyg with 
the VLBA and phased VLA. From their non-detection at 
15 GHz, when the VLA lightcurve suggested that an 8a 
source should have been seen, they inferred that the source 
was being resolved out, putting a lower limit on the source 
size of 1 mas. Together with light-travel time arguments ap- 
plied to the flares they detected, they constrained the source 
size to lie in the range 4-7 AU. While this is consistent with 
our upper limit of 1.3 mas, it might initially appear to be 
difficult to reconcile with the theoretical predictions of 0.3- 
0.6 mas. Either the models are out by a factor of 2-3 (quite 
possible given some of the inherent assumptions), or some 
unknown instrumental problem (possibly similar to the ap- 
parent suppression of the flux density we found when using 
the VLA during the flare) caused the source brightness to be 
reduced below their expected 8<r detection. In order to re- 
solve the jet, high-frequency, global VLBI would be required, 
to overcome the scattering (which scales as v~ 2 ) and achieve 
a resolution of ^ 1 mas. Should the jets not be detected 
with higher-resolution observations, it would present a seri- 
ous challenge to the jet interpretation of the radio emission 
from quiescent systems. 



5 SCATTERING 

The line of sight to the Cygnus reg ion is known to be highly 
scatter- broadened (|Fev et al.| [l989). Fluctuations in the elec- 
tron density in the interstellar medium (ISM) modify the re- 
fractive index of the plasma, inducing distortions in a wave- 
front travelling through the medium. This reduces the coher- 
ence of the electric field measured by two separated anten- 
nas, reducing the amplitude of the visibilities measured on 
long baselines, or, equival ently, giving ris e to angular broad- 
ening in the image plane. IRickettI i|l990t ) discusses other ef- 
fects of interstellar scattering, such as scintillation and tem- 
por al broaden i ng of p ulsed emission. 

IFev et aT] (|l989l ) measured the scattering properties 
along the lines of sight to several extragalactic calibra- 
tors in the Cygnus region, including our phase reference 
source, J 2025+3343 (i.e. B 2023+336) and our amplitude 
check source, J 2023+3153 (i.e. B 2021+317). The scatter- 
ing is clearly non-uniform, and the scattering properties can 
differ dramatically with small position offsets on the sky. 
Between our two calibrators, separated by only 1.9°, the 
scattering measure varies by a factor of five. The measured 
scattering size of J 2025+3343 at 4.99 GHz was 2.9 + 0.7 mas, 
scaling as u 1 , with 7 = -1.97 ± 0.29. At 8.4 GHz, we mea- 
sured a scattering disc size of 2.2 x 1.4 mas 2 , fitting in ei- 
ther the image or uv-plane. This is significantly larger than 
the previously-derived scattering size, although the size is 
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known to vary; iFev fc DesaJ (|200Ch claimed evidence for a 
50 per cent increase in scattering disc size at all frequencies, 
in observations taken three years after their original data. 

5.1 Source distance 

V404 Cyg l ies behind the Cygnus superbubble. 

lUvamker etail l|200ll ) found that the superbubble was 
not a single physical entity, but a projection of unrelated 
features at different distances. Nevertheless, the structures 
are all believed to lie between 0.5 and 2.5 kpc from us. 
While the line of sight to V404 Cyg does not pass through 
any of the known OB associations in the region, we can 
limit the location of the major scattering screen to be 
within the same range of distances. Assuming a single 
thin phase screen is responsible for the majority of the 
scattering, then according to the relative distance of the 
scattering screen and scattered source, we can predict the 
expected scattering size, 



(1 _ -fO^xgal, 



(3) 



where # ga i is the scattering size of the galactic source at 
distance d, # xg ai is the scattering size of the background ex- 
tragalactic source, and d s is the distance fr om the observer 
to the phase screen. As demonstrated by lLazio fc Cordesl 
ljl998t ). the observed scattering size of a galactic source de- 
pends on the relative distance of the source and the inter- 
vening phase scattering screen. The closer the source is to 
the scattering screen, the less scattered the source will be. 

The major axis of the calibrator scattering disc, 2.2 mas, 
together with the maximum major axis size of ~ f.3mas 
permitted by the fits to the source data in the image and uv 
planes shows that if the assumption of a thin phase screen 
is valid, and the scattering measure is similar along the lines 
of sight to the source and calibrator, then the lower limit to 
the source distance is d < 2.4d s . This gives limits of < 1.2 
and < 6.1 kpc respectively, for screen distances of 0.5 and 
2.5 kpc, consistent with the range o f 2.8-6.2 kpc range de- 
rived by Ijonker fc Nelemansl (|2004h . While this is not ter- 
ribly constraining, it shows that lower-frequency observa- 
tions which resolved the scattering disc of V404 Cyg could 
potentially provide an alternative, albeit model-dependent, 
estimate of the source distance. 



5.2 Scintillation 

Given the high scattering towards the source, we should 
determine whether the observed variability is intrinsic to 
the source, or whether it could be caused by interstel- 
lar scintillation. From the scattering measure of SMq = 
0.34 ± 0.12 m~ 20//3 kpc measur ed towards the phase refer- 
ence source bv lFev et al.l (|l989h . and assuming a scattering 
screen located at d p = 0.5 kpc, the n the standard scatter- 
ing equations (e.g. I Good man 1993) imply that we are in 
the strong scattering regime (for which the cutoff frequency 
is v < 100(rf s /d ) 5/17 (S'M/SMo) 6/17 GHz with these pa- 
rameters). Our observing bandwidth of 32 MHz is too large 
to detect diffractive scintillation (the decorrelation band- 
width is 1.6(d s /d )- 1 (5A//S'Afo)~ 6/5 MHz), and the source 
size, limited to > 5 fias by the inverse Compton limit of 
lb = 10 12 K, is also too large (diffractive scattering would 



require < 0.06(d s /rfo) _1 (5M/SM )" 3/5 //as). Thus we are 
in the regime of refractive scintillation. For intrinsic source 
sizes between 5 /ias and 1.3 mas (from the brightness tem- 
perature limit and the measured upper limit on the source 
size from the visibilities) , the timescale of refractive scintilla- 
tion is between 140 and 1000 h (scaling almost linearly with 
d s ), with a modulation index 0.02 < raR < 0.18 (scaling as 
~ (d s /do)- 1/6 (SM/SM ) 1/2 ). Thus the observed variabil- 
ity, an increase by a factor of 3 on a timescale < 1 h, must 
be intrinsic to the source. 



6 FLARING AND THE NATURE OF THE 
QUIESCENT STATE 

Despite being in a "quiescent" state, V404 Cyg exhibits 
sub-orbital v ariability in the o ptical (| Wagner et al.l I1992T ) 
and infrared (|Sanwal et a l. 1996) b ands, with some evidence 
for a 6 -h quasi-periodicity ( Casares et al.|[T993l ) . iHvnes et al.l 
(|2002h found that variable photo-ionization of the accretion 
disc by the X-ray source was responsible for the observed 
Ha flares, a conclusion reinforced by observations o f cor- 
related X-ray and Ha variability (|Hvnes et al.l [20041 ). The 
optical continuum f lares d i d not correlate as well with the 
X-ray flares, leading Irlvnesl (|2005l ) to suggest that the contin- 
uum variability could be due to synchrotron emission from 
a jet. The source is k nown to vary in the radio band (e.g. 
iHiellming et al.ll2000t ). and iMioduszewski et al.l (|2008T l pre- 
sented tentative evidence for a similar 6-h radio periodicity. 
The similar variability characteristics (timescales and flare 
amplitudes) seen at the different wavebands are suggestive 
of a common origin. 

While certainly intrinsic to the source, and clearly a 
common phenomenon in this system, the nature of the radio 
flaring is still to be determined. A thermal bremsstrahlung 
origin is ruled out, since it would produce too high an X- 
ray luminosity for the system. This suggests that the radio 
emission is non-thermal (certainly plausible given the rel- 
atively high brightness temperature), although a spectral 
index measurement during a flare is needed to determine 
whether such flares are optically thick or optically thin. Due 
to a paucity of simultaneous multiwavelength observations, 
any connection with the optical continuum flares remains 
purely speculative. We note, however, evidence for a correla- 
tion between the optical and radio bands is seen in the iden- 
tical power-law decay timescale measured in the two wa ve- 
bands during the 1989 outburst (|Han fc Hiellmmd[i"99"l ). 

The < 30 min rise time of the flare implies that the 
source of the flaring emission must b e of size < 4AU, ~ 2 
times the size of t h e bin ary orbit (|Shahbaz et al.l Il994l ). 
IMioduszewski et al.l (|2008l ) suggested three possibilities for 
the origin of such flaring; collimated jets, shocks in the cir- 
cumstellar medium, or diffuse emission as seen in SS 433 



(|Blundell et al.ll200ll) . The latter was on much larger scales 
(several hundred AU), and was interpreted by the authors as 
bremsstrahlung emission, for which reasons, we are inclined 
to favour the jet or shock interpretations. 

The fl at radio spectrum seen in the quiescent state of 
V404 Cyg (|Gallo et al.ll2005l ) implies optically-thin thermal 
or optically-thick synchrotron emission. Thermal emission is 
ruled out due to overpredicting the observed X-ray emission. 
Optically-thick synchrotron emission is typically interpreted 
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as a partially self-absorbed conical jet in X-ray binary sys- 
tems, which would make flaring in the jet, possibly due to 
internal shocks, a potential candidate for the observed vari- 
ability. Other circumstantial evidence for the existence of 
jets in the quiescent state is seen in the correlation between 
the X-ray and radio luminosities of hard-state black hole 
X-ray binaries at X-ray luminosities below L x ~ 0.01Z/Edd 
l|Gallo et alj|2003h . This correlation was later demonstrated 
to hold across the entire black hole mass spectrum (Merloni, 
Heinz & di Matteo 2003; Falcke, Kording & Markoff 2004), 
and simultaneous X-ray and radio observations of A0620-00 
ijGallo et al.l 1200(1 ) demonstrated that this correlation per- 
sists in quiescence, down to X-ray luminosities as low as 
10" 8 ' 5 L E dd. The resolved jets seen in GRS 1915+105 and 
Cygnus X-l at the high- luminosity end of the correlation 
then make it plausible that jets are also present at much 
lower luminosities. However, there are several caveats to this 
argument. A number of sources have recentl y been found 
with lower than expected radio luminosit i es dCorbel et al 



20041 ; IChatvll2006l ; ICadolle Bel et al.ll2007l ; iRodriguez et a 



2007), callin g into question the universality of the correlation 



( Galldl2007f ). Also, in the absence of direct imaging, nothing 
is known about the extent and degree of collimation of such 
quiescent jets, should they exist. The in ferred larger inner 
disc radii derived in the quiescent state (|McClintock et all 
2003) could potentially affect the collimation of the out- 
flow. While the SED of the quiescent source A0620-00 
has been well-fitte d with a maximally jet-dominated model 
l|Gallo et alj|2007f l. other models cannot be ruled out. Ex- 
tended, flat-spectrum radio emission not arising from a jet 
has been seen in the relativel y uncollimated expanding struc- 
tures observ ed in CI Cam ([Mioduszewski fc Rupen 120041 ') 
and RS Q ph (iRup en et alj|2007h and the equatorial emission 
of SS433 jBlundell et al.l 120011 ). Thus the true morphology 
of the emission in V404 Cyg cannot be definitively ascer- 
tained without high-resolution radio imaging. 

lHan fe Hiellmind (|l992h suggested that the short- 
timescale variability in the radio lightcurves of the 1989 
outburst of V404 Cyg could be interpreted as shocks 
at the ends of jets. Circumstellar shocks have been di- 
rectly imaged in the high-m ass X-ray binary CI Cam 
|Mioduszewski fc Rupenl l2004h . This system has a very 
dense stellar wind, which is thought to have smothered the 
jets ejected during the outburst, giving rise to an expanding 
shell of radio emission seen on VLBI scales. While the wind 
in V404 Cyg is likely to be less strong, the size constraint 
on the flares of V404 Cyg implies that they are occurring 
closer to the binary than the shock seen in CI Cam, where 
any wind from the donor star would be densest. With the 
available single-frequency, unresolved radio image, we can- 
not hope to pin down the location of the shocks. Ultimately, 
some explosive event is required to cause the radio flares, 
and whether they are shocks within or at the ends of jets 
cannot yet be determined. 

As the most luminous black hole X-ray binary in quies- 
cence, V404 Cyg is the only source in which such flaring ac- 
tivity has been detected in the quiescent state. The question 
of whether the flares are unique to this source, or are com- 
monplace in such systems, has implications for the nature of 
the accretion process at low luminosities. Owing to the lim- 
itations of current instrumentation, it has not been possible 
to observe other quiescent systems with sufficient time reso- 



lution in the radio band to detect such flares. Models fitting 
the SEDs of quiescent sources often assume a non-varying 
source, such that non-simultaneous observations in different 
band s could be used to build up the SED (e.g. iGallo et al.l 
If flaring emission is ubiquitous in these sources, the 
need for strict simultaneity will have to be more rigorously 
enforced. 

Observations of XTE J 1118+480 and GX 339-4 in their 
hard states (currently differentiated from the quiescent state 
only by a fairly arbitrary luminosity criterion) have shown 
evidence for fast variability in the opt ical and X-ray regimes 
|Motch et al.ll 19831 ; iHvnes et al.ll2003f ). In XTE J 1118+480, 
both the autocorrelations and the cross-correlation between 
the two wavebands were inconsistent with X-ray reprocess- 
ing and more indicat ive of optical synchrotron variability 
|Kanbach et all 1200 if ), which Malzac, Merloni & Fabian 
(2004) were able to explain with a jet-dominated model. 
The timescales are shorter in this source than seen in V404 
Cyg however, and there were no high-time resolution radio 
data for comparison (although we note that we would not 
expect to see radio variability on timescales shorter than a 
few minutes, since any variations occurring at the base of a 
jet would be smoothed out further downstream where the 
source became optically thin in the radio regime). A bet- 
ter comparison is the quiescent black hole in the centre of 
our Galaxy, Sgr A*, accreting at ~ 10 _8 LEdd- It shows ra- 
dio and infrared flaring activity (although any connection 
between the two bands is as yet unclear), and the radio 
flaring has been explained as adiabatic expansion of a self- 
absorbed transient pop ulation of relativistic electrons (e.g 
lYusef-Zadeh et alJl2006tY 

There is thus significant evidence that flaring activity 
in or at the end of a jet is a feature of the quiescent state 
of accreting black hole systems. Whether there is a truly 
steady underlying jet as a ssumed by standard jet models 
(IBlandford fc Konig]|[l979l ) or whether the emission is com- 
posed of multiple overlapping flares, or some more uncol- 
limated structure, remains to be determined. More sensi- 
tive instruments such as the Expanded VLA (EVLA), e- 
MERLIN and the upgraded VLBA are necessary to probe 
these short-timescale flares and determine the nature of the 
quiescent emission. 



7 CONCLUSIONS 

We have detected the quiescent black hole X-ray binary sys- 
tem V404 Cyg with the High Sensitivity Array at 8.4 GHz. 
The unresolved source has a brightness temperature of > 
10 7 K, and the emission is inferred to be non-thermal syn- 
chrotron radiation. Our detection provides the most accu- 
rate source position to date, which can serve as an initial 
point for future proper motion studies. We have also con- 
strained the length of the jet in this system to be < 1.3 mas, 
consistent with previous measurements and theoretical ex- 
pectations. 

A small flare was detected with both the VLA and the 
HSA, with a rise timescale of 30min, in which the source 
flux density rose by a factor of 3. This flare is certainly in- 
trinsic to the source, and cannot be caused by interstellar 
scintillation. This is further evidence that source flux den- 
sities in the quiescent state are not stable, mandating the 
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use of strictly simultaneous observations when constructing 
and modelling SEDs to determine the contribution of differ- 
ent components (such as jets) to the observed emission in 
different wavebands. 

The detection of an unresolved quiescent black hole X- 
ray binary system opens the way for further high-resolution 
observations. If the inferred jets indeed exist, and are flat- 
spectrum, as has been proposed, then higher-frequency ob- 
servations would be well-placed to resolve the jets, if the 
theoretical predictions for the jet length are correct. Fur- 
thermore, the unresolved source at 8.4 GHz is an ideal can- 
didate for a parallactic distance determination, free from the 
uncertainties associated with the source structure seen when 
such systems are observed with VLBI during outburst. 
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